The Annals of Applied Probability 
2008, Vol. IS, No. 4, 1569-1587 
DOI: 10.1214/07-AAP496 

© Institute of Mathematical Statistics, 2008 

DISORDERED PINNING MODELS AND COPOLYMERS: BEYOND 

ANNEALED BOUNDS 

By Fabio Lucio Toninelli 1 

Ecole Normale Superieure de Lyon and CNRS 

We consider a general model of a disordered copolymer with ad- 
sorption. This includes, as particular cases, a generalization of the 
copolymer at a selective interface introduced by Garel et al. [Euro- 
phys. Lett. 8 (1989) 9-13], pinning and wetting models in various 
dimensions, and the Poland-Scheraga model of DNA denaturation. 
We prove a new variational upper bound for the free energy via an 
estimation of noninteger moments of the partition function. As an 
application, we show that for strong disorder the quenched critical 
point differs from the annealed one, for example, if the disorder dis- 
tribution is Gaussian. In particular, for pinning models with loop 
exponent < a < 1/2 this implies the existence of a transition from 
weak to strong disorder. For the copolymer model, under a (restric- 
tive) condition on the law of the underlying renewal, we show that 
the critical point coincides with the one predicted via renormaliza- 
tion group arguments in the theoretical physics literature. A stronger 
result holds for a "reduced wetting model" introduced by Bodineau 
and Giacomin [J. Statist. Phys. 117 (2004) 801-818]: without restric- 
tions on the law of the underlying renewal, the critical point coincides 
with the corresponding renormalization group prediction. 

1. Introduction. We consider a rather general class of directed poly- 
mers interacting with a one-dimensional defect through a quenched disor- 
dered potential. This includes various models motivated by (bio)-physics: 
among others, wetting models in (1 + 1) dimensions [12, 14], pinning of 
(1 + d) -dimensional directed polymers on columnar defects [27], copolymers 
at selective interfaces [16, 25] and the Poland-Scheraga (PS) model of DNA 
denaturation [11, 24]. For further motivations and references, we refer to 
[17], Chapter 1. One of the interesting aspects of these models is that they 
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present a nontrivial localization-delocalization phase transition due to an 
energy-entropy competition. 

Mathematically, the model is defined in terms of a renewal sequence whose 
inter-arrival law has a power-like tail with exponent a. + 1 > 1. The model is 
exactly solvable in absence of disorder, and it turns out that the transition 
can be of any given order, from first to infinite, according to the value of 
a. This is therefore an ideal testing ground for physical arguments (Harris 
criterion, renormalization group computations) and predictions concerning 
the effect of disorder on the critical exponents and on the location of the 
critical point. 

The comprehension of this model has witnessed remarkable progress on 
the mathematical side, as proved by the recent book [17]. In particular it 
has been shown that for wetting, pinning or PS models, an arbitrary amount 
of disorder modifies the free-energy critical exponent if a > 1/2 [18], that 
is, disorder is relevant in this case, in agreement with the predictions of the 
so-called Harris criterion [23]. On the other hand, for < a < 1/2 it has been 
proven recently [2, 29] that if disorder is weak enough the free-energy crit- 
ical exponent coincides with that of the homogeneous (i.e., nondisordered) 
model, and the (quenched) critical point coincides with the annealed one: 
disorder is irrelevant (again, in agreement with the Harris criterion). These 
results about "irrelevance" of disorder for < a < 1/2 have been later re- 
fined and complemented in [21] with results about correlation-length critical 
exponents. The marginal case a = 1/2 is strongly debated in the theoreti- 
cal physics literature: Ref. [14] claims that quenched and annealed critical 
points coincide for disorder weak enough, while [12] concludes the opposite 
and gives a precise prediction for their difference. See [2] and [29] for rigorous 
results in the marginal case, which however do not solve the controversy. 

Here we attack two major open problems: 

• Do quenched and annealed critical points coincide for strong disorder! 
The Harris criterion, which is based on the analysis of the stability of 
the homogeneous model to the addition of weak randomness, makes no 
prediction about this point for pinning models with a. < 1/2 or for the 
copolymer with any a. [Here and in the following, we say for brevity 
"pinning models" but we actually include wetting and PS models, besides 
pinning of (1 + d) -dimensional directed polymers on columnar defects. 
Mathematically all these are variants of the same model, cf. beginning of 
Section 3.1.] 

• For the copolymer model, it is known that the critical point is bounded 
above by the annealed one and below by an a-dependent expression found 
by non-rigorous renormalization group arguments [25]. Are either of these 
two bounds optimal? 
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In this work we prove a new upper bound on the free energy of the model 
which in some cases is sufficient to answer these two questions. In particular, 
consequences of our bound include the following: 

1. Both for pinning and copolymer models with, say, Gaussian randomness, 
for large disorder quenched and annealed critical points differ. We would 
like to emphasize that, especially for the "marginal case" of the wetting 
model with loop exponent a = 1/2, this question was subject to dispute 
even very recently [12, 15]. 

2. We identify the strong-disorder behavior of the critical point both for 
Gaussian pinning models and for a "reduced" wetting model introduced 
by Bodineau and Giacomin [5] (for the "reduced model," the same result 
was proven recently by Bolthausen, Caravenna and de Tiliere [6]; their 
proof is however based on a very different method). 

3. For the copolymer model we prove that, as soon as a homogeneous depin- 
ning term is present in the Hamiltonian, quenched and annealed critical 
points differ for every strength of the disorder; in particular, the much- 
studied critical slope at the origin (cf. Section 3.3) is in this case strictly 
smaller than 1. 

4. Finally, again for the copolymer model we prove that, if the law of the 
underlying renewal sequence satisfies a certain explicit condition [cf. equa- 
tion (3.43) below; in particular, the condition requires the renewal to be 
transient], the critical point predicted by nonrigorous renormalization 
group arguments is indeed the correct one. This is however believed not 
to be the case in general, that is, if (3.43) does not hold. 

Our basic idea is to estimate noninteger moments EZ" 7 of the partition 
function with l/(l+a)<7<l. The reason why we cannot go down to 
7 < 1/(1 + a) is not just technical and will become clear soon. The method 
of fractional moments has been already applied successfully to the study of 
other quenched disordered models: we would like to mention in particular 
(a) the work [1] by Aizenman and Molchanov, where bounds on small mo- 
ments (of order less than 1) of the resolvent kernel of random Schrodinger 
operators are employed to prove the occurrence of Anderson localization for 
strong disorder or extreme energies; (b) Ref. [8] by Buffet, Patrick and Pule 
who compute exactly the free energy of a directed polymer in random en- 
vironment on a regular tree via an estimation of EZ 7 with 1 < 7 < 2; and 
(c) Ref. [13] where Derrida and Evans, again via an estimation of EZ 7 with 
1 < 7 < 2, improve previously known estimates on the critical temperature 
of directed polymers in random environment on finite-dimensional lattices. 

We would like to conclude this introduction with two remarks. First of all, 
for pinning models with < a < 1/2 and Gaussian randomness our results, 
together with those of [2] or [29], imply that there is a nontrivial transition 
from a weak- to a strong-disorder regime, see Remark 3.3 below for a precise 
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statement. Secondly, Theorem 3.6 together with the numerical simulations 
of Ref. [9] strongly indicates that the critical point of the copolymer depends 
not only on a but also on the details of the inter- arrival law of the renewal, 
a possibly nonintuitive fact. 

2. The model and the main result. Let r := {to, t\, . . .} be a renewal se- 
quence on the integers, with inter-arrival law K (•): tq = and {ji — Tj_i}j>i 
is a sequence of IID random variables taking values in NU {+oo}, with law 
P(ri = n) =: K(n). We assume that for n € N 

(2.1) K(n) = ^ 

with < a < oo and L(-) a function varying slowly at infinity, that is, a pos- 
itive function such that lim-r-^x, L{rx) / L(x) = 1 for every r > 0. In general 
E re6N ^(n) = l-P(n = +oo)<l. 

A very popular example in the literature is the one where K(-) is the 
law of the first return to zero of the one-dimensional simple random walk 
{S k }k>o, that is, K{n) = K SRW (n) := P(inf{£; > : S k = 0} = 2n\S = 0) (in 
which case a = 1/2, Xm>i K{n) = 1 and L(-) ~ const). 

In order to make a quick connection with the (bio)-physical systems 
presented in the introduction, let us mention that in the case of (1 + d)- 
dimensional pinning models one takes a = d/2 — 1 if d > 2 and a = 1/2 if 
d = l, in (1 + l)-dimensional wetting models and copolymers at a selective 
interface usually a = 1/2, while a ~ 1.15 in the case of the PS model in three 
dimensions [24]. Actually, for the copolymer model and the pinning model 
with d= 1 one usually makes the specific choice K(-) = K SRW {■). 

We consider random copolymers with adsorption. The system has size 
N € N, and it is characterized by the parameters /3>0, A>0, /i€M and 
h > 0. Moreover, we let u := {w„} ne N and uj := {uj n } n ^ be the realizations 
of two IID sequences of random variables. We assume u to be independent 
of uj. The joint law of (u;,tD) is denoted by P. 

The partition function of the model is 



Z N,u;,Z '■- E 



(2.2) 



{n6i) 

2 J MN€T} 

where In(t) := |r n {1, . . . ,N}\ = J2n=i l{nGr} 5 while the free energy is 
(2-3) F((3,h,\,h) ^lim^logZ^-. 
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We assume as usual the existence of all exponential moments of uo\: 
(2.4) M{u) :=E(e UUJ1 ) < oo 

for every u£l, and similarly for M(u) := E(exp(ucDi)), and we set by con- 
vention Ewf = Kujf = 1 and EcDi = (this can always be achieved via a 
redefinition of f3, A and h). For later convenience (cf. Section 3.2), we do 
not assume in general that the u; n 's are centered, although this can always 
be obtained via a trivial shift of h. Under these assumptions on P, it is well 
known that the limit in (2.3) exists almost surely and in L 1 (P), and that it 
is almost surely independent of (u,u>). Another well-established fact is that 
F(f3, h, A, h) > 0, which is an immediate consequence of 

J3u! N +h 

(2-5) Z NjU ~ > - K(N) 

and (2.1). For a proof and a discussion of these facts, see for instance [17], 
Chapters 1 and 4. 

It is actually one of the main questions in this context to decide when the 
free energy vanishes and when it is positive. The reasonjs the following: One 
usually defines the localized region as C := {(/3,h,X,h) : F((3,h, \,h) > 0} 
and the delocalized region as T> := {((3,h, X,h) : F(f3,h, X,h) = 0}. Then it 
is well known that, if (f3,h,X,h) G £, the ratio ~(In(t))/N converges 
for N — ► oo to a positive limit, almost-surely independent of (uj,u) (where 
Ejy w ~ denotes the disorder-dependent Gibbs average) , while the same limit 
is zero in the interior of T>. This explains the names given to the two regions. 
We refer to [19] and [17], Chapter 8, for more precise and more refined 
estimates on In(t) i n the interior of T>, to [4], [20] and [17], Chapter 7, for 
path properties of r in C and finally to [28] for estimates on In{t) at the 
boundary between V and C. 

A very cheap way of showing that the system is delocalized for given 
values of (/?, h, A, h) is via the annealed upper bound on the free energy: 

F(p,h,X,h) 

<F ann (p,h,X,h) 

(2.6) 

:= lim — logE 

TV^oo N 



e I N (r)(h+logM(l3)) 



/ 1 + e (-2hA+logW(-2A))(T i -r J -_i) x 

x II [ 2 JVer} 

which is simply Jensen's inequality: ElogZ < logEZ. If we define the delo- 
calized region of the annealed model as V ann := {(/?, h, A, h) : F ann ((3, h, A, h) = 
0}, one has then obviously V ann C V. 
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The point of this work is to provide an improved upper bound on F which 
is enough to conclude that D ann ^ T>. In some cases, we will even be able to 
identify sharply the boundary between the two regions. 

Going beyond annealing has appeared so far to be a difficult task. A 
natural idea is to try Morita-type bounds [26], that is, constrained annealing. 
In other words, for every function p(uj,u>) such that Ep(a>, u) = 0, it is easily 
seen that 

(2.7) F((3, h, X,h) < lim i-logE(e^)z ~). 

JV — >oo iv 

This can indeed improve the upper bound (2.6) on F if p is suitably chosen 
but, as shown in [10], it cannot improve the estimate T> ann C V, as long as 
p is a sum of local functions ( "finite-order Morita approximation" ) : 

N 

p(uj,uj) = ^2po(6 ti (uj,lj)), 

n=l 

with po(u>,lu) a bounded function depending only on a finite number of 
(ui,uji), and 9 being the shift operator (8 n (uj,uj)) m = (uj n+m ,uj n+rn ). We will 
show that the estimation of noninteger moments of the partition function 
allows to bypass this difficulty. 

In order to formulate our results we need a few auxiliary quantities. For 
< 7 < 1 let 

(2.8) c( 7 ) := £ (K{n)Y< > c(l) = 1 - P(n = +oo), 

neN 

with strict inequality if 7 < 1. We remark that c(-) is decreasing and that 
0(7) <ooifl/(l + a)<7<l. Also, observe that 

C (i/(i+a))=x;— — 

neN n 

can be finite or infinite according to the behavior at infinity of L(-). 
Define also, for n 6 N and 7 such that 0(7) < 00, 

{K(n)y< 



(2.9) K^(n 



c{i) 



so that XmeN-^7( n ) = 1 by the definition of 0(7). It is important to realize 
that Kry{-) is still of the form (2.1), just with a replaced by (1 + a)j — 1 > 
and L(-) by L(-) 7 /c(7) (which is still slowly varying at infinity). 

Finally we need the following definitions: for a, b, v S R, k € N and < 

7<1. 



(2.10) f 7 (a,b;k):=E 
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(2.11) 



G~j(v, a, b) 



:= lim — log EL 

N->co N 7 



In(t) 

e ui N (r) -Q / 7 (o,6;r i -r i _i)l {JVeT} 



where P 7 is the law of the (recurrent) renewal with inter-arrival law K- f {-). 
The limit exists by super additivity and is nonnegative. 
Our main result is the following: 

Theorem 2.1. Let 1/(1 + a) < 7 < 1 be such that 0(7) < 00. Then 

(2.12) F(B,h,X,h) < -GJlogch) + /17 + logM(/3 7 ), -2A, -2AM. 

7 

Some interesting consequences of this result are worked out in Section 3. 

PROOF of Theorem 2.1. We have the elementary: 

Lemma 2.2 ([22], Chapter 2.1). Let 0<7<1. IfneN and ax > 0, . . . , 
a n > 0, then 

(2.13) (ai + --- + o n ) 7 <a 7 + --- + a 7 . 
We need also the identity 



N 



(2.14) 



E E 



=10=io<*l<— <i*=JV \j=l 
£ 

><n 

i=i 



JJ ^(ij- - j exp ^9 £j Wi . + M^J 



l + exp[-2A£^ = , i+1 (£ n + /0] 



which is just a way of rewriting the expectation in (2.2) as a sum over all 
possible configurations of r n {1, . . . , N}. As a consequence of Lemma 2.2, 

N / i \ 



E E n*e 

Le=lO=i <h<--<i e =N \j°=l 



(2.15) 



x exp 



f^E 



* ^ 1 + eX p[-2A E^_ 1+ i (2„ + /»)] x 
x II 
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Performing the disorder average, one sees that the right-hand side of (2.15) 
equals 



E E (ri^-^i)V 

i=10=in<ii<— <it=N \7=1 / 



£[log M (7/3)+ft-7+log c(7)] 



^=10=i <ii<— <fy=N \j= 

(2.16) 



Therefore, 



(2.17) 

< G 7 (logc( 7 ) + /i 7 + logM(/? 7 ), -2A, -2Xh). 



limsup — ]ogE(Z N ~y 

N->oc 



On the other hand for 7 > we have via Jensen's inequality: 

(2.18) F(P,h,X,h)= lim i-ElogZ ~ < limsup -i- log E(Z ~p 

which concludes the proof of Theorem 2.1. □ 

3. Applications to disordered pinning models and copolymers. 

3.1. Random pinning model. In this section we assume that A = 0, Ku\ = 
and we call for simplicity the free energy F(/3, h). Since A = the random 
variables u) n do not play any role, and for simplicity we write Zn,w instead 
of Z N u ~ for the partition function. 

The model thus obtained is then the one considered for instance in [2, 
12, 14, 18, 29]. According to the law P and, especially, to the value of a, 
the model is interpreted in the physics literature as a pinning, or wetting, 
or Poland-Scheraga model in different spatial dimensions. 

We observe that F((3, •) is nondecreasing and we denote as usual by h c ([3) 
the (quenched) critical point of the pinning model: 

(3.1) hjp) :=w£{heR:F(p,h) > 0}, 

while the function j3 1— > h c {0) will be referred to as the critical curve. When 
(3 = (homogeneous pinning model) it is a standard fact that h c (0) = 
— logP(ri < +00) (cf., for instance, [17], Chapter 2): F(0,h) is positive 
for h > h c (0) and zero otherwise (for the detailed behavior of F(0,h) for 
h \ h c (0), cf. [17], Theorem 2.1). It is also known (see [3] and [17], Chapter 
5.2) that h c {(3) < h c (0) for every (3>0. 
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The annealed bound (2.6) applied to this case shows that 
h c ([5) > K nn ((3) := -logM(/3) - logP(n < +00) 

(3.2) 

= - logM(p) -logc(l). 
On the other hand, since / 7 (0,0; k) = 1, Theorem 2.1 implies immediately 



Theorem 3.1. For every f3 > one has 

(3.3) hc(J3)>h c (P):= sup --log[M( 7/ 3)c( 7 )]. 

l/(l+a)<7<l 7 

Of course, /i c (-) depends on K(-) because c( 7 ) does. 
The important point is that the bound provided by Theorem 3.1 is in 
various cases strictly better than the annealed one. For instance: 

Corollary 3.2. Assume that logM(/3) af3 p for some p > 1 and 

a > [where A(x) x ~°° B(x) means lim x ^ OQ A(x)/B(x) = 1/. 
Then, there exists (3q < oo suc/i i/iaf /or enery (3 > (3o 

(3.4) ^09) >/ir n (/3). 

This applies for instance to the centered Gaussian case lo\ ~ JV(0, 1) where 
logM(/3)=/? 2 /2. 

Proof of Corollary 3.2. Choose a value of 7 € (1/(1 + a),l). One 
has for (3 — > 00: 

(3.5) - - log[M( 7 /3)c( 7 )] ~ -wf- 1 ^ 

7 

while h^ nn (P) ~ —af3 p , and the statement is obvious from the conditions 
p > 1 and 7 < 1. □ 

One can easily extract from Theorem 3.1 a sufficient condition which 
guarantees that h c (/3) > h^ nn {j3) for a given /3 > 0, that is, 



(3.6) a 7 (--log[M( 7 /3)c( 7 )] 



<0. 

7=1 



For instance, in the case of Gaussian disorder and recurrent renewal [i.e., 
c(l) = 1], this condition can be re-expressed in the simple form: 

(3.7) ^>-]>>(n)logA». 

1 n=l 
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In the recent work [15] it is claimed, on the basis of numerical simulations, 
that in the case of the (1 + l)-dimensional wetting model [which corresponds 
to a = 1/2, L(-)~ const and c(l) = 1/2] with to n taking only two values 
with equal probability, quenched and annealed critical points coincide even 
for strong disorder. Theorem 3.1 does not disprove this assertion because for 
symmetric two- valued w n 's it turns out that equation (3.3) does not improve 
the annealed bound, but in our opinion it makes the scenario suggested by 
[15] rather unlikely. 



Remark 3.3. Corollary 3.2 is particularly interesting when < a < 1/2 
and disorder is Gaussian. Indeed, together with the results of [2] or [29], 
it implies that there exist < 0\ < 0% < oo such that h c (0) coincides with 
the annealed critical point for < 0\ and differs from it for > 02- I n this 
sense, one can say that a transition from a weak disorder regime to a strong 
disorder regime occurs. 



3.1.1. Strong- disorder asymptotics of the critical point. It is clear that, 
under the assumption of Corollary 3.2 on M{0), for very large the supre- 
mum in (3.3) is realized by some 7 very close to 1/(1 + a). This observa- 
tion, together with a generalization of ideas from [5], allows to identify the 
strong-disorder asymptotic behavior of the critical point. For instance, in 
the Gaussian case one has: 

(3.8) h c {0) 



2(1 + a) 



and it is easy to obtain an analogous statement in the case \ogM{0) 
a(3P. 



Proof of equation (3.8). It is immediate to see from (3.3) that 
(3.9) Hminf^#> 1 



02 - 2(1 + a)' 

As for the opposite bound, it is based on a straightforward generalization 
of the rare stretch strategy of [5] (for this reason, we just sketch the main 
steps of the proof). Let £ € N, assume that N is an integer multiple of £ and 
divide {1,2, . . . ,N} into blocks I k := {(k - 1)1 + 1, (k - 1)1 + 2, . . . , k£}, with 
k = 1,2, . . . , (N/£). Given q > and the disorder realization uj, let 



(3.10) T w := ll < j < (N/£) : £ u n > £q \ U {N/£} 
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One obtains a lower bound on the partition function as follows: 

Z Nu} > E[e£«=i (/3w " +/l)1 {"^>; 

(3.11) 

h c r VfcG^rnA = 0V^I 1J ,k (JV/i)]. 

In other words, we have constrained r to visit each point in the last block 
and in each of the blocks where the empirical average of the w n 's is larger 
than q, and to skip all the others. We can now take the logarithm, divide by 
./V and let N — > oo at £ fixed in (3.11) to get a lower bound on the free energy. 
We do not detail this step, since an essentially identical computation appears 
in the proofs of [5], Proposition 3.1, [18], Theorem 3.1 and [17], Theorem 
6.5. The net result is that for every e > 



(3.12) F((3,h)>p(£) 
where 

(3.13) p{£) : 



(3q + h + \ogK{\) - (1 + a + e)^- + o/(l) 



and oi(l) is a quantity which vanishes for £^> oo. The term log-ftT(l) is due 
to P(/iCr)=K(l/. 

From (3.12) and the definition of the critical point one deduces that 

(3.14) h c (P) <-(3 q + (l + a + £ )^-- log Jf(l). 

Optimizing over q and using the arbitrariness of e > 0, 

(3-15) h c (l3)<-—^—-\ogK(l) 

2(1 + a) 

which, together with (3.9), proves equation (3.8). □ 

3.1.2. About the size of Z^,u in the delocalized phase. For the homoge- 
neous model, it is known [17], Theorem 2.2, that if h < h c (0) then 

(3.16) Z N = E[e hI ^ T h {Ner} ] N ^°° CK(N) 

where C > depends on h and on P(ti < oo). It is natural to ask whether 
a similar statement holds for the disordered model inside the delocalized 
phase. In this respect, the ideas developed in this work allow to go much 
beyond the statement of Theorem 3.1 that the infinite- volume free energy 
vanishes for h < h c ((3), and prove the following: if h < h c {(5) there exists 
1/(1 + a) <7<1 and a constant C := C(h,K(-)) < oo such that 

(3.17) W(Z Nu , > uK(N)) < Cu-i 
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for every u > 0. The upper bound (3.17) on the partition function should be 
read together with the lower bound (2.5): both are of order K(N), just like 
the estimate (3.16) which holds for the pure model. 

PROOF of equation (3.17). Since h < h c ((3), there exists 1/(1 + a) < 
7 < 1 such that 

A := - log[M(/?7)c(7)] - hy > 0. 

Then, it follows from the proof of Theorem 2.1 [cf. in particular equations 
(2.15) and (2.16) taken for A = 0] and from (3.16) that 

(3.18) E(Z N)W y < E 7 [e-^M A l {Arer} ] < CK 7 (N) = JL-K{Nf, 

for some C := C(h,K(-)). Equation (3.17) is then an immediate consequence 
of Markov's inequality. □ 

Let us remark also that one can extract from the proof of (3.17) the fol- 
lowing almost sure statement: if h < — (I/7) log[JVf(/?7)c(7)] for some 1/(1 + 
ck) < 7 < 1, then 

(3.19) JNjuN^oo 
P(da->)-almost surely, for every [i < 1 — 1/(7(1 + a)). 

3.2. " Reduced" wetting model. This model, introduced in [5], Section 
4 as a toy version of the copolymer model discussed in the next section, 
is obtained from (2.2) putting A = h = and assuming that ¥(u>\ = 1) = 
p = 1 — P(cji = 0), for some p E (0, 1). Moreover, one assumes here that the 
renewal is transient, c(l) < 1, otherwise the model is uninteresting in that no 
phase transition occurs. As we will see, the actual value of c(l) is irrelevant 
in the strong-disorder regime we are going to consider, as long as it is strictly 
smaller than 1. Changing conventions with respect to the previous section, 
we denote the free energy as F{(3,p) in this case. 

It is known that for every (5 > there exists p c (P) E (0, 1) such that the 
free energy F(f3,p) is zero for p < p c ((3) and positive for p > p c (j3). The main 
question is to compute (if it exists) the limit 

(3.20) m c := - lim -logp c (/3). 

p— >oo p 

It was proven in [5] that for a = 1/2 

(3.21) 2/3 < -limsup^logp c (/3) < - liminf - logp c (j9) < 1. 

/3-»oo P P 
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This result can be easily generalized to any a > and in this case the lower 
bound in (3.21) is replaced by 1/(1 + a). 
Here we prove the following: 

Theorem 3.4. The limit in (3.20) exists and equals 1/(1 + a). 

Recently a proof of Theorem 3.4, based on entirely different ideas, was 
given in [6]. While our approach is much simpler, the method of Bolthausen, 
Caravenna and de Tiliere is more natural from a renormalization group point 
of view. 



Proof of Theorem 3.4. We start by remarking that if we apply The- 
orem 2.1 to the reduced model we obtain immediately 

(3.22) -liminf-logp c (/3)<l J 

j3— >oo p 

but not the sharper statement of Theorem 3.4. To go beyond (3.22), somehow 
we have to use the information that, if p < p c (P), the sites where u) n = 1 are 
very sparse for (3 large (their density p is indeed exponentially small in /3), 
and that between two such sites r has typically just a finite number of points, 
since it is a transient renewal under the law P. 
Let 1/(1 + a) < 7 < 1 and C > 1 such that 

(3-23) ^ < 1. 

The renewal r with inter-arrival law K{-) being transient, that is, c(l) < 1, 
it is known (cf., for instance, [17], Theorem A. 4) that 

(3.24) P(n G r) n ~°° _!_#(„), 

(l-c(l)) 2 

and therefore there exists C := C"(7, C) < oo such that 

(3.25) P(ner)<C"^ 
for every n £ N. 

We let now := {1 < n < N : uj n = 1} U {N} and we decompose the par- 
tition function according to the configuration of A(t) := r n y u . If |A(r)| = 
1)> we write A(t) = {ai,a2, ■ • ■ ,ae} with m < aj+i and, by convention, 
we set ao := 0. Then, 

n e 

(3.26) Z NJU <Y, E e<i " II P (°i " a J-i G r ) 

^=1 ^C^: j = l 

\A\=£,a e =N 
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N 



(3.27) <]T Y, e^ +lo ^ e Y[ 

1=1 ACy u : j=l 
\A\=£,a e =N 



(/3+iogC'K TT K ( a i ~ Q j-i) 
C 



(3.28) 



N e 
<e /3+logC'^ e (/3+logC')JJ j=1 ^. 

1=1 0=a <a 1 <---<a t =N 

x -Q -g(flj -Qj-l) 



where in the first inequality we used the fact that 

P(a, Gr,rnXn {a. t _i + 1, . . . ,04 - 1} = 0|oi_i £ r) 

(3.29) 

< P(oi - a,_i € r), 

in the second one we used (3.25) and the third one is obvious since we have 
just added extra positive terms to the sum. On the other hand the right- 
hand side of (3.28), apart from the global factor exp(/3 + logC") which is 
anyway irrelevant for the computation of the infinite-volume limit of the 
free energy, is just the partition function of the model where @ is replaced 
by (/3 + logC) and K(n) by K{n)/C for every raeN [cf. equation (2.14) 
taken for A = 0]. Following step by step the proof of Theorem 2.1, one finds 
therefore that 

lim ^-ElogZ N ^ 

N-+00 Is 

(3.30) <limsup-^logE(^) 7 

< lim iiogi^W^t*)^)^ ^?-?))!! , 
which vanishes whenever 

(3.31) logfre^W) + (1 -p)} +log(^) < 0. 
Thanks to the choice (3.23), it follows immediately that 

(3.32) -liminf-logp c (/?)< 7 . 

p— *oo p 

By the arbitrariness of (1 >) 7 > 1/(1 + a), and since we already know 
that 

-limsup^logp c (/3) > 1/(1 + a), 
we obtain the statement of the theorem. □ 
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3.3. Copolymer model. In this section we set /3 = h = and with yet 
another abuse of notation we call the free energy F(X,h). This is just (a 
generalization of) the copolymer model considered for instance in [4, 5, 7]. 

Remark 3.5. The results which follow can be easily generalized to the 
case h ^ 0. This is particularly evident for h < 0. Indeed, it is well known 
[and immediate to check from (2.14)] that the model with h< is equivalent 
to the one with h = and K(n) replaced by K(n) exp(— \h\) for every n £ N 
[so that 0(7) becomes c(7)exp(— 7| 



We observe that in view of A > the free energy is nonincreasing with 
respect to h and we denote by h c (X) the quenched critical point: 

(3.33) h c (X) :=sup{h>0:F{X,h)>0}, 

while A 1— > h c (X) is the critical line. 
It is convenient to define for < 7 < 1 

(3.34) /,(7) (A):= JL log M(-2A 7 ). 

Note that h { c\x) is increasing by Jensen's inequality and that 

(3.35) l im MLW =7 
y J A\0 A 

thanks to Kiof = 1, Eu3i = 0. For 7 = 1, this is just the annealed critical 
line: F ann (0, 0, A, h) = if and only if h > hc(X), as it is immediate to 
realize from (2.6). On the other hand, /ic 1+ (•) is sometimes referred to 
as the "Monthus line." This line was proposed as the true critical line of 
the copolymer model in the theoretical physics literature, on the basis of 
renormalization group arguments [25]. 

In general one knows that for every A > 0, 

(3.36) ^(V (i+«)) (V) < h c (X) < hP (A) . 

The lower bound in (3.36) was proven in [5] in the case where K(-) = 
K SRW '(■), the law of the first return to zero of the one-dimensional sim- 
ple random walk (cf. Section 2), but the bound can be proven to hold for 
every K{-) of the form (2.1), see [17], Chapter 6. Numerical simulations 
(supported by solid probabilistic estimates) performed in the particular case 
K(-) = K (■) strongly indicate that neither of the two bounds in (3.36) 
is optimal [9]. 

For K (•) = K SRW (•), it is also known that the limit slope lim^o h c (X)/X 
exists [7] [and is therefore bounded between 2/3 and 1 in view of (3.36) and 
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a = 1/2], that it coincides with the slope of the critical curve of a continuous 
copolymer model with Brownian disorder [7] and that it is to a large extent 
independent of P [19]. 

Here we will show that (say, in the Gaussian case) the annealed upper 
bound h c (X) < h^\\) can be improved for large A whatever K(-) is [within 
the class (2.1)]. If the renewal is transient, Z)neN-^( n ) < li then with no 
assumptions on the disorder distribution annealing can be improved for ev- 
ery A > and, in particular, the slope at the origin turns out to be strictly 
smaller than 1. Finally, if K(-) satisfies condition (3.43) below, the lower 
bound in (3.36) is the optimal one for every A. This condition is not sat- 
isfied for K(-) = K w {-), in agreement with the simulations mentioned 
above. By the way, our results strongly indicate that the critical curve de- 
pends in general not only on the tail behavior of K(-) (say, on the exponent 
a), as one might be tempted to guess on the basis of belief in universality, 
but also on the details of the slowly- varying function L(-). Note that this 
is not the case for the annealed curve, which depends only on the disorder 
distribution P, nor for the Monthus line which depends only on P and a. 

Observe first of all that, thanks to Lemma 2.2, 

1 + Mby+logM(ja)) 

(3.37) / 7 (a, b; k) < 2 1 ~T— . 

Therefore, 

G 7 (v, a, b) 

(3.38) < lim ^logE 7 e %M("+(i-7)i°g2) 

n + g(b7+logA/( 7 a))(r J -T j _ 1 )x 
11 !{iVGr} 



X 

Our main result for the copolymer model is the following: 

Theorem 3.6. Define 

(3.39) 7c := inf{0 < 7 < 1 : log 0(7) < 7 log 2} G [1/(1 + a), 1). 
Then, one has: 

1. For every 7 > 7 C there exists C(j) < 00 such that for every A > 

(3.40) £ c (A)<fcW(A) + £M. 
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2. J/ 7 C = 1/(1 + a) and 

a 

(3.41) < log c(l/(l + a)) + log 2 < log 2 

1 + a 

there exists C < oo suc/i £/iai /or every A > 

(3.42) Uc(A) <fc(V(i+«))(A) + ^. 

A 

3. If 7c = 1/(1 + a) and 

(3.43) iog C (i/(i + a)) + __ log2<0, 

1 + a 

iften /or every A > 

(3.44) ne(A) = n.( 1/(1+a)) (A). 

In view of Remark 3.5 above, condition (3.43) is realized for instance if 
c(l/(l + q)) < oo and we add a homogeneous depinning term proportional to 
— \h\, with \h\ sufficiently large. In any case, we emphasize that a necessary 
(but not sufficient) condition for (3.43) to hold is that r be transient under P. 

Remark 3.7. Note that in the case L(-) ~ const, for example, if K (•) = 
K w (-), one has 7 C > 1/(1 + a) since 0(7) — > 00 for 7 \ 1/(1 + a). In 
particular, from the explicit expression of K SRW (-) one finds numerically 
7 C ~ 0.83 in the simple random walk case. It is interesting to note that the 
numerical results of [9], Section 5, are compatible with h c (X) = he {X) with 
m somewhere between 0.8 and 0.84. Understanding whether this is more 
than just a coincidence requires further numerical simulations of the type 
[9], performed also for other choices of K{-). 

We have also the analogue of Corollary 3.2: 

Corollary 3.8. Assume that logM(— A) ~ aA p for A — > +00, for some 
a > and p > 1 . Then, there exists Ao < 00 such that for A > Ao 

(3.45) h c (X)<h^(X). 
This is obvious from (3.40) since j c < 1. 

Finally, as we mentioned, in the transient case (and without assumptions 
on P) the improvement on annealing can be pushed down to A = 0: 

Corollary 3.9. Assume that c(l) = P(ti < +00) < 1. Then, for every 
A > 

(3.46) h c (X)<hP(X) 
where 

(3.47) 7 := inf{ 7 < 1 : log 0(7) + (1 - 7) log2 < 0} € [1/(1 + a), 1). 
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As a consequence of equation (3.35) we have that if P(ti < +00) < 1 (or 
if h < 0, see Remark 3.5) the slope of the critical curve at the origin (if it 
exists) is strictly smaller than 1. 

Proof of Theorem 3.6. One has from (3.38) 
G 7 (logc( 7 ),-2A,-2A[7 l W(A)+ e ]) 



(3.48) 



< lim — logE^ 

_ N->oo N • 



3 / JV (T)(log C ( 7 )+(l-7)log2) 

/1 +e -2A7£(T J -T 3 _ 1 ) 

k n 



l{Ner} 



Let us consider first case (1). To this purpose, let 7 > 7 C and choose e = 
(7/(7 A) in (3.48). It is clear that, for every 5 > 0, it is possible to choose 
C = C{5) sufficiently large so that 

(3-49) ±±±. 



< e 



lo g 2+<5 



for every k € N. Taking for instance 
(3.50) 



J = J(7) = 7 Jog2-logc(7) >0 



(so that C depends on 7) one finds therefore 



(3.51) 



G 7 logc(7),-2A,-2A 



C(7) 



7 A 



< lim llogE 7 [e-^)^V er} ] = 0. 

[The positivity of ^(7) follows from 7 > 7 C and from the definition of 7c .] By 
Theorem 2.1 one has therefore (3.40). Equation (3.42) is proven analogously. 

Finally, if condition (3.43) is verified then choosing 7 = 1/(1 + a) and 
e = in equation (3.48) one has 

G m+a) ( log c( T ± T ) , -2A, -2\h^ 1+a » (X) 



(3.52) 



1 + a 



1 



< lim ^-logE 1/(1+Q) [e^W( lo s c ( 1 /(i+«))+«/(i+«)io g 2) 1{7vg7 



and the right-hand side is zero since it is the free energy of a homogeneous 
pinning model with nonpositive pinning strength, see assumption (3.43). 

This shows that h c (X) < hc^ 1+ (A), and the opposite bound is already in 
(3.36). □ 
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Proof of Corollary 3.9. Just take equation (3.48) for 7 = 7 and 
e = 0. □ 
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